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There were two principal objectives of the cooperative agreement between NASA and the 
University of Denver. The first goal was to modify the design of the ER-2 condensation 
nuclei counter (CNC) so that the effective lower detection limit would be improved at 
high altitudes. This improvement was sought because, in the instrument used prior to 
1993, diffusion losses prevented the smallest detectable particles from reaching the 
detection volume of the instrument during operation at low pressure. Therefore, in spite 
of the sensor’s ability to detect particles as small as 0.008 microns in diameter, many of 
these particles were lost in transport to the sensing region and were not counted. 

Most of the particles emitted by aircraft are smaller than 0. 1 pm in diameter. At the start 
date of this work, May 1990, continuous sizing techniques available on the ER-2 were 
only capable of detecting particles larger than 0.17 micron. Thus, the second objective of 
this work was to evaluate candidate sizing techniques in an effort to gain additional 
information concerning the size of particles emitted by aircraft. 

Accomplishments during this Cooperative Agreement 

Redesign of the ER-2 CNC . The design of the ER-2 CNC was evaluated under this grant 
and revisions were made in the design which would permit a larger fraction of small 
particles to reach the detection volume of the instrument during low pressure operation. 
These design changes were implemented in the construction of the ER-2 CNC II which 
was used in the 1993 Stratospheric Photochemistry, Aerosol and Dynamics Expedition 
and in ASHOE/MAESA flights which are currently under way. Figure 1 shows the design 
changes incorporated into the ER-2 CNC as a result of this work. The redesigned 
instrument, the ER-2 CNC II, was constructed with support from another grant. The 
response of the current, redesigned CNC is less dependent on pressure than that of the 
previous instrument. 

Evaluation of sizing techniques. A number of techniques were evaluated for determining 
the sizes of particles smaller than 0.1 micron in diameter. These techniques included a 
differential mobility analyzer (DMA) and a diffusion battery (DB). In addition, the 
possibility of extending the lower size range of the optical particle counter was also 
considered. The object was to provide accurate size distributions of particles emitted from 
aircraft flying in the stratosphere. 

Evaluation of the DMA technique 

Technical reports generated in the study of the DMA are attached in Appendices 1 and 2. 
A summary of the results is given here. 

The DMA technique promises high resolution measurements of aerosol size distributions. 
The technique consists of exposing the aerosol to a bipolar ion atmosphere which will 



neutralize most of the particles but leave a certain number with a single charge, fewer with 
a double charge, etc. After this conditioning, the aerosol is introduced into the DMA 
where particles are exposed to known flow and electric fields. These fields are crossed, so 
that only particles within a narrow range of electrical mobilities can reach the exit region. 
Thus, ignoring the effects of multiple charging, particles within a very narrow size range— 
virtually monodisperse-exit the DMA. These particles are counted by a CNC. By 
varying the voltage and/or the flow, a spectrum particle concentration vs. particle mobility 
may be measured. Accounting for multiple charging effects, this spectrum can be 
converted to a size distribution for particles larger than the CNC detection limit (about 
0.008 jam diameter). Thus, this technique hold the promise of high-resolution size 
distribution measurements of particles in the size range likely to be produced by jet 
engines. 

The analysis of the DMA alternative in Appendices 1 and 2 shows that there are two 
fundamental issues which remain unresolved for a stratospheric DMA. The first issue 
involves the effects of Brownian diffusion on the transfer function of the DMA. 
Discrepancies between the predicted and actual transfer function for very small particles— 
those most affected by diffusion— exist in the literature. These discrepancies suggests that 
the effects of diffusion are not precisely predictable or experimentally controllable. Since 
diffusion is increasingly important as pressure decreases, the response of a stratospheric 
DMA may not be fully predictable. The magnitude of this effect has not been ascertained. 

A second problem associated with the operation of the DMA at stratospheric pressures 
involves the charging of aerosol at low pressure (Appendix 2). The recovery of a size 
distribution from a DMA depends critically on the distribution of charges on aerosol 
particles resulting from the ion preconditioning process. Although ion properties and 
aerosol chargeing have been studied at atmospheric pressure, there is little literature on 
these processes at low pressure. It may be possible to use the preexisting ambient charge 
distribution and avoid the the artificial preconditioning problem. At the current time, the 
knowledge of stratospheric ion concentrations, ion-attachment coefficients at low 
pressure, and the charge distribution on particles found in aircraft exhaust is insufficient to 
evaluate the effectiveness of a stratospheric DMA. 

A third problem involving a stratospheric DMA, one not discussed in the appendix, is 
associated with the absolute concentration of particles available in the stratosphere. Only 
a small fraction of the particles upstream of the DMA are charged in the preconditioning 
process. Since the concentration of stratospheric particles is quite small, obtaining a 
reasonably accurate size distribution measurement in the stratosphere would require 
extremely long sampling times to obtain a statistically significant signal. The aerosol being 
measured would have to remain relatively constant over the long measuring time. Even in 
the plume of an aircraft, the large spatial variations in particle concentrations would make 
complete particle size measurements with a DMA impractical. 

Given the uncertainties in diffusion losses within a DMA at low pressure, the unknown 
distribution of charges on stratospheric particles, and the low concentrations of these 



particles that imply very long sampling times, using a DMA for measuring stratospheric 
aerosol size distributions appears to be impractical. Some of these uncertainties may be 
reolved with future laboratory measurements and modelling studies; however, the problem 
of long sampling times will probably require fundamental changes in DMA design and may 
remain unsolved. 

Evaluation of diffusion battery technique 

Diffusion batteries have been used for decades to evaluate the size distribution of small 
particles in the atmosphere. In this class of instruments, particles are separated by passing 
the aerosol through a “battery” of screens or tube arrays. Particles with high diffusional 
mobilities (i.e., small particles) diffuse to the walls of the battery and are removed from the 
airstream. The remaining particles are counted by a CNC. 

There are several problems associated with determining particle size distributions from 
diffusion batterie measurements. The main difficulty lies in the nature of the response 
function and the difficulty in recovering accurate particle size distributions from the 
response function. These difficulties are discussed in more detail in Appendix 3. Briefly, 
the diffusion battery does not provide a sharp discrimination of particle size, rather a very 
gradual one. The measured parameters are strongly convoluted with the response 
function and must be deconvolved, or inverted, to obtain the original distribution. 

Because this convolution is of the form of a Fredholm integral equation of the first kind, 
there are an infinity of possible solutions (particle size distributions) that could produce 
the measured response. Various techniques are availabe that constrain the equation and 
produce the “best” solution, while another reports all physically possible solutions, giving 
a wide range of particle size distributions. When measurement errors are added to the 
inversion, the recovery of a particle size distributions becomes very difficult. Because of 
these problems, DMAs are diffusion batteries are rarely used now in cases where a DMA 
can be applied. 

Because of the fundamental problems associated with inverting diffusion battery data and 
recovering a size distribution, it is unlikely the diffusion batteries can be used to evaluate 
the details of the size distribution of small particles in the stratosphere. However, the 
technique shows some limited promise in providing rough estimates of some important 
parameters— such as mean particle size— in the diameter range from 0.008 to 0.1 pm. 

Redesign of the PMS aerosol spectrometer 

In response to requests made by Dr. Guy Ferry of NASA Ames Research Center, Particle 
Measuring Systems Inc. (PMS, Boulder, CO) offered to build an optical particle counter 
capable of sizing particles as small as 0.05 pm in diameter. Since the instrument heats 
particles before measuring them and thus reduces their size, a sulfuric acid particle having 
a diameter of 0.05 pm at the laser would have a diameter approaching 0.07 pm under 
ambient temperatures and water vapor mixing ratios. This instrument, the focused- cavity 
aerosol spectrometer, or FCAS, represents a considerable advance on the previous model, 



the active-scattering aerosol spectrometer. NASA has purchased the upgraded optical 
particle counter. The FCAS flew successfully on the AASE-II and SPADE rfiissions 
before being replaced with a model with an improved laser and new pulse analysis 
hardware. The FCAS instruments have collected data that indicate the presence of a mode 
of small particles when the CNC counts high concentrations of particles, although this 
mode is only partially resolved by the spectrometer measurements. 


Summary 

The ER-2 CNC has been modified to reduce the diffusive losses of particles within the 
instrument. These changes have been successful in improving the counting efficiency of 
small particles at low pressures. 

Two techniques for measuring the size distributions of particles with diameters < 0.17 jam 
have been evaluated. Both of these methods, the DMA and the diffusion battery, have 
fundamental problems that limit their usefulness for stratospheric applications. We cannot 
recommend either for this application. Newly developed, alternative methods for 
measuring small particles include inertial separation with a low-loss critical orifice and 
thin-plate impactor device. This technique is now used to collect particles in the multi- 
sample aerosol collector housed in the ER-2 CNC-II, and shows some promise for particle 
size measurement when coupled with a CNC as a counting device. 

The modified FCAS can determine the size distribution of particles with ambient diameters 
as small as about 0.07 pm. Data from this instrument indicate the presence of a nuclei 
mode when the CNC-II indicates high concentrations of particles, but cannot resolve 
important parameters of the distribution. 


Bibliography. 

No publications resulted from these engineering and feasibility studies. 

Figure Captions: 

Figure 1. Assembly A shows the design of the ER-2 CNC. In this instrument, the sample 
was drawn through the sample flow meter, Al, up the injector in the vertical saturator, 

A2, and into the condenser. In the revised design. Assembly B, the flow meter and 
injector are combined which leads to a shorter length from the inlet to the saturator. Thus 
fewer particles are lost by diffusion. 

Figure 2. Modelled response function of a tube-array diffusion battery under stratospheric 
conditions. 
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Figure 2. Modelled response function of a tube-array diffusion battery under stratospheric 
conditions. 
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STRATOSPHERIC DMA 


NOMENCLATURE 

Dimensional variables: 
t — time 

r — radial coordinate 
z = axial coordinate 
ft = flow velocity at r,z 
E = electric field at r,z at time t 
V = analyzer voltage at time t 

Dimensional parameters: 

Z p = particle electric mobility 

D = particle diffusivity 

Ri = radius of central rod 

R 2 — inner radius of outer cylinder 

Dh = 2 (R 2 -R\) = hydraulic diameter 

L - axial distance between aerosol entrance and exit slits 

A z a = axial width of aerosol entrance slit 

A Zs = axial width of aerosol exit slit 

Q a = aerosol entrance flow 

Q s = aerosol exit flow 

Q c = sheath entrance flow 

Q m = sheath exit flow 

Qt = Qa+<2c = Qs+Qm = total flow 

U = Qt/^-Ri 2 ) = mean flow velocity 

Vq = analyzer voltage at f=0 

t = analyzer voltage scan time constant 

Zp* = (<? ra +Qc)-log(Ri/R 2 )/47rTK= centroid of steady-state transfer function at V 
Z p0 . = Z p -(V=V 0 ) 

Dimensionless variables: 

1 = t/r 
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u={r/R 2 ) 2 
z = z/ -Dh 


Dimensionless parameters: 

7 = (i^i/ R2) 2 
/C = LR2/{R2 7 -Rl 2 ) 

«a = Az a /D 
ACs = A^s/I 

ft = ( Qs+ 5a)/ ( 5m + 5c) 

S=(5s-5a)/(5s+5a) 

G = geometric flow factor 

I) = AirLD/ ( <2m+ 5c) 

^ = G-I) 

(7*2 — (j2(^p=^p*) = G'D * 

Re - pUDh/fi = Reynolds number- 
5= L/rU 
Zp — Zpf Zpo* 


INTRODUCTION 

It is proposed to develop a new particle sizing system to be operated aboard the 
ER-2 in the stratosphere. One of the possible measurement systems consists of a 
differential mobility analyzer (DMA) operated upstream of a condensation nucleus counter 
(CNC). The purpose of the DMA is to extract aerosol within a known narrow size range 
and present it to the CNC for concentration measurement. The response of the DMA is 
well understood in its standard steady-state mode (Knutson and Whitby, 1975a). 

However, for this project it proposed to operate in the scanning mode described by Wang 
and Flagen (1990). Use of the scanning mode may present several problems depending on 
the rate of scan. Some of these problems may be confounded by particle Brownian 
diffusion which is significant for the proposed application. Aerosol charging is also a major 
area of uncertainty. 

In the discussion below, the DMA scanning mode is first analyzed neglecting 
diffusion. Diffusion effects are discussed in the following section. Aerosol charging is 
discussed in a separate paper. 
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SCANNING MODE AND CHARACTERIS TIC TIMES 

Strictly laminar flow may be assumed within the DMA classification section itself. 

In general, laminar flow may be assumed from some indeterminate point upstream of the 
DMA all the way to the CNC viewing volume. There may be some mixing at certain 
__ points in the transport lines, such as just inside the DMA aerosol exit slit, which will affect 
the following analysis. But for the most part, particles may be assumed to follow flow 
streamlines in the absence of diffusion and the DMA electric field. 

In analyzing the transit of a particle through the proposed measurement system, 
there are several key locations along its trajectory from the sample point to the detector. 
The first is where it enters the sampling system at the tip of the external sample probe, 
designated as p. Well mixed flow is assumed from p to b, the point where laminar flow 
begins upstream of the DMA. The DMA aerosol entrance and exit slits are designated as a 
and 5 , respectively, the CNC detection volume as d. Let r m represent the average particle 
transit time through the well mixed zone from p to b. Let r a , r x and r s represent the 
average particle transit times from b to a, a to 5, and 5 to d, respectively. Let tp, £b, k, k, 
and td represent the times at which a particular particle passes the respective points. An 
inherent assumption for the well mixed zone is that all particle transit times are the same, 
It. t b -fp=rm. But for the laminar flow regions the transit times for a particular particle 
depend on which streamline it is on. Because of slow flow near walls, the variation of 
transit times through the laminar zones upstream and downstream of the DMA is on the 
order of the mean transit times, r a and r s . The variation of transit time through the DMA 
itself is significantly less than r x . Let r represent the period over which the voltage of the 
scanning DMA changes by a significant amount, and n the period over which the ambient 
aerosol size distribution changes by a significant amount. The CNC sample interval is 
designated as ra- 

Even when operated in the steady-state mode, the size range extracted by the DMA 
has a finite width determined by the ratio of aerosol to sheath flows. When operated m the 
scanning mode there are several factors acting to broaden the particle size range presented 
to the CNC at any given instant. 

Particles arriving at the CNC at a given instant have exited the DMA at point s 
over a range of earlier times. If r s «r does not hold, then the mean particle size exiting 
the DMA has changed significantly over this time range resulting in a broadening of the 
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size range presented to the CNC at that instant. 

If t x «t , then the DMA voltage may be treated as constant over the period of a 
single transit of the DMA classification section and there is no broadening of the extraction 
range for this part of the particle trajectory. However, as the scan rate is increased (r 
decreased) the DMA transfer function is increasingly warped from its basic triangular 
shape and becomes broadened. (See results of SCAN simulation program.) 

If rd<<r does not hold, then again the particle size exiting the DMA has changed 
significantly over the sample interval resulting in broadening of the effective average 
extraction range over the interval. 

For normal stratospheric sampling Tf may be assumed to be large compared to the 
sampling system transit time. This condition may be strained somewhat for the High 
Speed Research Program when flying through a fresh jet plume. If the laminar flow zone 
transit time is not small compared to Tf then particles arriving at the CNC at a given 
instant have been sampled from the ambient over a range of times during which the 
ambient size distribution has changed significantly. 

In general, if r and Tf are large compared to T a , t x> t s and r<i then the 
measurements can be treated as steady-state. If T>>Td does not hold, then the response of 
the measurement system is still relatively easy to analyze in terms of an average size 
extraction range over the sample period. If t is not large compared to T a , r x or t s then the 
problem becomes much more difficult. Particle diffusion makes theoretical analysis of the 
effect of transit time variation extremely difficult. 


DMA TRANSFER FUNCTION IN SCANNING MODE 

Consider a particle traversing the classifier from the aerosol entrance slit to the 
aerosol exit slit. It starts at radial position u= 1 on a streamline designated by ft n and 
ends at radial position u= 7 on a streamline designated by ft u t- ftn and ftut are 
proportional to stream function. ft n is 0 on the center streamline of ft, -1 on the 
downstream edge of the entrance slit and +1 on the upstream edge. ftut is 0 on the center 
streamline of ft, -1 on the downstream edge of the exit slit and +1 on the upstream edge. 
In the steady-state mode, the dimensionless mobility ft of a particle traversing from 



(1 , Oin) to (7, Qout) does not depend on the flow velocity profile. This is not the case in the 
scanning mode. 


Unless the flow profile can be accurately determined by numerical modeling, the 
DMA transfer function must be calibrated experimentally. The dimensionless flow velocity 
profile b/ Uis a function of position (cj,z) and depends on the dimensionless parameters Re , 
A 5 , 7, *a and %. Calibration is required over the entire range of these parameters 
anticipated to be encountered under ER— 2 flight conditions. For a given DMA, 7, k, 
and Ks are fixed by geometry. /?, 6 are determined by the relative values of the DMA flows. 
If the flow split is kept constant then /?, 5 remain constant. Re is likely to vary during 
flight due to the large pressure range covered. A constant mass flow control strategy would 
be required to keep Re constant. This is difficult in view of issues discussed below. 

If the functional time dependence of the scanning DMA voltage is of the form 
V= Vo-exp(t/r) where ris positive or negative, then all particles traversing the classifier 
from (l,<2in) to (7,§ 0 ut) follow the same trajectory. In this case, DMA transfer can be 
characterized by the function IF= W(Qi m Qout,S) where Sis the dimensionless scan rate. If 
7is not of exponential form then the transfer characteristic depends on particle mobility 
(size) as well, making experimental calibration much more arduous. 

Assuming an exponential form for V, the sensitivity of the transfer characterization 
function W to the flow profile, hence, to Re, /?, 5 , 7, k, k & and /Cs, can be estimated by 
calculating IF for several trial profiles. The basic profile used consists of pure axial flow for 
the bulk of the region between u=l and u= 7 except for pure radial flow regions at the 
aerosol entrance and exit slits connecting streamlines in the entering and exiting aerosol 
flows to the corresponding streamlines in the axial flow region. Slug and Poiseuille axial 
profiles were used matched with uniform radial profiles in both cases. A uniform radial 
electric field was assumed. 

Particle trajectories were calculated for the above set of flow and electric fields. A 
particle entering on streamline at u=lj2—1 follows that streamline radially (at 
increased velocity due to radial electric field) to radial position a/i n (0in) w T here it enters the 
axial flow region. It moves axially due to the flow and radially due to the electric field 
across fluid streamlines until it reaches streamline Q 0 ut at radial position w ou t(0out)- II 
then follows that streamline radially (at increased velocity) to radial position 
where it exits the classification section. and w 0 ut are determined by 
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Q(^in) = (H?in)/5(l-^)/2(l+/?) 

Q(w 0 ut) = 1 ” (l+Oout)/5(i + ^)/2(l+/5) 

where 

$M = ihy- / 1 s*(w' )d «' 

1 UJ 

is the fraction of total flow between radial position w and 1, and ih.{u)=v-J U \s the 
dimensionless axial velocity profile. 


V7is the solution to the equation 

U{n u z (w) 

XW-u ) 




&0Ut 


do; . 


The dimensionless mobility Z p of the particle traversing the classifier from (1, Qin) to 
(7,§out) and located at radial position Wi n >a)>w 0 ut in the axial flow region at time t is given 

by 

wr = w+ [ir?]-l l - eVT ■ 

If U, tin, tout and h are the times the particle is at radial positions a 2 , Win, Wout and u h 

respectively, then these are related by 

(« 2 -lin)/r = W-ui 1_. : 

W-bJ\n tV-U in [ 

p(^ouf^in)/ 7 ’ W-»ov\ 

C W-U in 

e (fl“tout)/^ _ -rr r-i 

v W- Wont VY-Uovt 

Derivations of these equations are given in the handwritten notes. 


h-t 


in 


^aj 

P+rr 

L ^ J 

[1-7 

(/}( 1 +<5)1 

^1“ t out 

„ «s. 

1+ ft 

* St 

V. 


The DMASCAN Fortran library is intended for use as an extension of the DMA 
Fortran library which is internally documented. DMASCAN contains the following 
functions with correspondence to the above notation as noted. 

FDENS(T,P) = p(T,p) _ 

ZSCAN(Q$IN,QSOT,STIM,LPRO,IND,TSCAN) = [Z p e >' ]{Qin,Qout,S,„{t r h)/T) 

EXPLIN(A.B) = root of e x +Bx-A 
WSCAN(WIN,WOT,GAM,STIM) = H'{win,Wout,7i5) 

DSCANL(W.VF) = [ih{u)-{W-u)/ 

WFLOW(F.G) = w(<?) 

VFLOW(W.G) = M“<7) ■ 
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The ZSCAN function returns the parameter ~Z p e'^ r which can be used to calculate the 
mobility of a particle which traverses the indicated path and arrives at the exit slit (u=U\) 
at time t\. LPRO is an integer flag indicating the axial flow profile: 0=slug, 2=Poiseuille. 
LPRO is passed to WFLOW and VFLOW through COMMON block FLOPRO. IND is an 
integer index referring to DMA configuration parameters stored in COMMON blocks 
defined by a call to DMASET in the DMA library. These COMMON blocks supply the 
required values of 7, 0 and 8 to ZSCAN. ZSCAN uses EXP LIN to solve the 

transcendental equations for the classifier transit times. The total dimensionless classifier 
transit time (ti-t 2 )/r is returned as an output argument of ZSCAN. WSCAN finds W 
through a search routine which calls TSIMP, a Simpson numerical integration routine in 
the DMA library, for which DSCANL calculates the integrand. WFLOW supersedes the 
DMA library version for applications requiring LPRO— 0. For the DMASCAN version to 
supersede the DMA version in linking, the DMASCAN library must be listed before the 
DMA library. FDENS is not used by any of the other DMASCAN functions but is used 
elsewhere for calculating Re. FDENS complements FVISC and FPATH in the DMA 
library which calculate air viscosity and mean free path. 


SCANF and SCAN are applications programs which utilize the DMASCAN library. 


Both programs are currently set up to give characteristics of the TSI Model 3071 DMA 
operated in the scanning mode at one atmosphere and 300K with 2 and 20 1pm aerosol and 
sheath flows. Both routines ask the user for LPRO and S. SCANF also requires Q\ n and 
O out and returns w in , w out , W, zJ'l T and (t r h)/T. SCAN returns Z p e k ' T and (t r U)/r 


for Qi n =0 and § O ut=0 and deviations of these parameters for an array of Q\ n and Qout 
values specified by the user. The values are printed on the screen and into the file 
SCAN.DAT. Files of the form PnSx.DAT are sample results of SCAN output. 


The actual DMA transfer function Q as defined in Wang and Flagen (1990) is 
obtained by some sort of integration of Z p e^l T over Q\ n and Qout- This has not yet been 
derived. 

If mixing of the flow streamlines does not occur between the classifier aerosol exit 
slit and the CNC detection volume, then a rigorous analysis of the effective instantaneous 
DMA transfer function as seen at the CNC would require taking into account the varying 
transit time t^-ts between the DMA and the CNC for each exit streamline Qout* 
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PARTICLE DIFFUSION 

Brownian diffusion of particles is an important effect for the proposed measurements 
of fine and ultrafine aerosols in the stratosphere. The particle diffusion coefficient is 
independent of pressure in the continuum regime and inversely proportional to pressure in 
the kinetic (free-molecule) regime. With a mean free path in air (at 23 C) of 0.067 fim at 
one atmosphere increasing to 1.36 fin i at 50 mbar, the aerosol of interest is mostly in the 
kinetic and transition regimes. Consequently, particle diffusion is of increasing importance 
in the higher altitude range of the ER-2. 

There are two major effects of particle diffusion on the operation of the DMA on the 
ER_2. The first is particle loss to the walls of sampling and transport lines. For laminar 
flow these losses can be predicted using the theory of Gormley and Kennedy (1949). The 
second effect is a broadening of the response function of the DMA. This is due to diffusion 
within the actual classification section of the DMA and in the transport tube from the 
DMA aerosol exit slit to the detection volume of the CNC. This second effect is much 
more difficult to analyze and is also coupled with the loss effect. 

Cross-stream diffusion causes particles to travel randomly between adjacent 
streamlines in the classifier and transport tubes. For a particle entering a transport tube 
on a given streamline, this leads to a probability distribution of transit times and exit 
streamlines. The author is not aware of any analytic or even numeric or Monte Carlo 
solution to this problem. Particle loss to the walls of the transport tube should have some 
limiting effect on the breadth of this probability distribution. Streamwise diffusion 
contributes also contributes to the distribution of transit times. 

A rigorous analysis of diffusion effects in the classifier is even more difficult. 
However, Kousaka et al (1985, 1986) and Stolzenburg (1988) have analyzed the problem of 
diffusion in the steady-state DMA. Kousaka et al have done a dimensional analysis of the 
problem and show results of an approximate numerical model and experiments. 

Stolzenburg has constructed an approximate analytic model which is perhaps more 
convenient for estimating diffusion effects in the stratospheric DMA. 


The non-diffusing DMA transfer function Q{Z P ) is triangular in shape (Knutson and 
Whitby, 1975a). As diffusion increases in importance, 0 becomes smaller and broader and 
more Gaussian in shape. The analytic model of Stolzenburg shows that [JQ(Z p )dZ p ] 
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remains constant at /3(l+5) while the mean (centroid) and standard deviation of fi (Z p ) are 
given by l+<r* 2 and [/? 2 (l+« 2 )/6+<r* 2 (l+<r‘ 2 )] 1 ' 2 , respectively. The dimensionless diffusion 
parameter o' is given by o' 7 =G- D\ D‘ is the dimensionless diffusivity of the particle 
corresponding to Z p = 1 and is the same as the parameter D/E extracted in the dimensional 
analysis of Kousaka et al. (1985). G is a geometric flow factor evaluated as an integral 
along the particle trajectory. For the TSI Model 3071 with 0=1/10, 6=0 and the same 
simplified flow field as described above, £7=4.1104 for Poiseuille flow and £7=3.5110 for 

slug flow. 

The steady-state analysis of diffusion in the DMA can be used as a first 
approximation for the diffusion effects in the scanning DMA. As the relative scan rate 5 
increases the error in this approximation also increases. The formulas above for the mean 
and standard deviation of fl (Zp) can be used to estimate the degree of the diffusion effect. 
The full form of the steady-state transfer function for diffusing aerosols can obtained and 
employed in data, reduction using the DMA library. 

The experimental data obtained by Stolzenburg (1988) to verify his diffusion theory 
indicate a discrepancy which may be significant for the stratospheric DMA. Figure 4.8 in 
Stolzenburg (1988) plots the relative error of the theory in predicting the mean mobility 
extracted by the second DMA in a tandem DMA (TDMA) experiment. Based on the . 
experimental measurements, the theory shows an apparent growth of the aerosol m going 
from DMAl to DMA2 when, in fact, there is no physical explanation for such growth. 

Other investigators have observed this same effect to a lesser degree for larger particle 
sizes. As seen from the top plot in figure 4.8, the effect increases as particle size decreases. 
The two experiments at different flow rates indicate that the effect also depends on 
something other than particle diameter. The lower plot shows a good correlation for both 
experiments with a measure of the relative broadening of the DMA transfer function due to 
diffusion. If this correlation is not merely coincidental it may indicate an error in the 
diffusion model which is proportional to the relative effect of diffusion. Depending on 
operational parameters, the relative diffusion effect in the stratospheric DMA is likely to 
extend beyond the range indicated in the lower plot of figure 4.8. Thus, substantial errors 
could occur in sizing the aerosol if the diffusion theory is in error as suggested. 

The experimental data of Stolzenburg indicate that the theory accurately predicts 
the width (standard deviation) of the transfer function. The discrepancy described above, 
may be due to an unpredicted skewing of the transfer function. There is some ambiguity in 
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the derivation of the form and use of G in the equations. The programs SCANG, GVAR 
and FVCO'R represent attempts to estimate the magnitude of the effect of these 
ambiguities on the theory. They indicate that the magnitude is insufficient to explain the 
discrepancy in the TDMA experiment, ihe effect of space charge in the TDM A 
experiment was also analyzed and also appears to be insufficient to explain the discrepancy. 
Hence, the discrepancy cannot be explained and could be significant for the stratospheric 

DMA. 


The DMASET subroutine in the DMA library evaluates G for the particle 
trajectory from Qi n =0 to Qout = 0 for Poiseuille flow. (Program GCON, not in the library, 
back calculates Gfrom variables in the DMA COMMON blocks.) The SCANG program, 
similar to SCAN, calculates G for an array of Q\ n and Qout values for either slug or 
Poiseuille flow. Contained within the SCAN.FOR file are functions FGEO and GINT. 
FGEO calculates Gfor specified Q- m and 5 0 ut values, flow profile flag (LPRO) and DMA 
configuration index (IND). FGEO calls GINT which supersedes the DMA library version 
for applications requiring LPRO=0. 
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AEROSOL CHARGING FOR STRATOSPHERIC 
TYn?rFRTvMTT AT, MOBILITY AN AL\ ZER 


INTRODUCTION 

A major area of uncertainty in reducing data from a differential mobility analyzer 
fnMA’l nnerated in the stratosphere is the question of the electric charge distribution on 

“! 0 Se £5^“^ 

toward this goal. A list of key investigators, one for each group, is. Fuchs, Hoppel, 

Marlow, Adachi, Pui, Wiedensohler, and Vienna (Porstendorfer et al.) 

The most promising theory to date for determining ion-aerosol attachment 
coefficients appears to be that of Fuchs. There seems to be general agreement that Fuchs 
theory works well in the continuum regime and well into the transition regime. This is t 
basis for most of the theoretical work of the groups mentioned above. Hoppel and Marlow 
have made the greatest efforts to modify Fuchs’ theory to obtain better agreement wit 
experimental data in the near free-molecule regime. The best data for making comparisons 
appears to be that of Adachi, Kousaka and Okuyama (1985) and of Pui, Fruin and 
McMurry (1988) though the latter is plotted in a much more convenient manner. Adachi 
compared his results to Fuchs’ original theory while Pm compared 5® p 
original theory and to a theory of Marlow. In a subsequent comment < 1 1 
Marlow and McFarlane (198S) showed improved agreement with Pui s data using a 
modified theory. In another comment, Hoppel and Frick (1989) showed good 
their theory with Pui’s data and state (from personal communication with Marlow) that 
there is an error in Marlow’s comment paper. Thus, Fuchs) theory as modified by Hoppel 
appears to be the best theory to date. What follows is a brief ° ver ^ ° f h 

work Parameter notation for the most part follows that used by Hoppel. Much or 
Hoppel ’s work relies heavily on previous work of others which is * ly h ^ erenced 1 

his papers. For the most part these references are not specifically mentioned here. 

There are five major papers for Hoppel’s work. Hoppel (1977) provides the 
foundation for calculating the ion-aerosol attachment coefficients. Hoppel liy&oj 
introduces (in his work) the timeffiependent ion-aerosol balance equations and shows some 
example solutions. Hoppel and Frick (19S6) give probably the most complete overall 
description of the aerosol charging problem. This paper reiterates most of the content o 
Hoppel (1977) and includes tables of attachment coefficients for typical tropospheric ions at 
nonnal [empeTture and pressure. It also demonstrates several solutions to the population 

balance equations including possible aerosol-aerosol q tera f 10n P he o n °TdnVv H a?affist the 
Frick (1989) show probably the best test of the attachment coefficient theory against the 
experimental uUrafine unipolar aerosol charging data of Pui, Fruin and McMurry (1988). 
Some information about the sensitivity of calculated attachment coefficients I to 110 
properties is provided. Hoppel and Frick (1990) treat the t.meKiependent .on-aerosol 
balance equations in depth, particularly with regard to the Cha ° 

distributions from neutralizers. Several solutions are shown for different neutralize 
configurations. The appendix provides a brief review of the problems involved in 
calculating attachment coefficients particularly with regard to^ the uncertainty ™ 
determining ion properties. Note is made of an improvement in .the ! fons 

calculation, and tables are included of attachment coefficients for typical neutralizer 10 . 
at normal temperature and pressure. The reader is cautioned that there are man. 
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typographical errors m the equations of these papers of which only a fraction have been 
discovered and noted. 


POPULATION BALANCE EQUATIONS 

The first (last) step in defining the aerosol charge distribution calculation is the 
system of time-dependent ion-aerosol population balance equations Each equation equates 
the time derivative of a particular species to a sum of terms describing the birth and death 

processes for that species. 

There is one equation each for positive and negative ions. Each ion equation has a 
generation/birth term (usually from radiation), an ion-ion recombination/death term, and 
death terms associated with elch possible charge state of an aerosol particle .There may 
also be death terms associated with losses to containment surfaces.. If more than one 
distinct ion species for a given polarity is considered (e.j. .negative : ion: i and free electrons), 
then there is one equation for each species and corresponding additional birth/death terms 
in all equations. Unless otherwise noted, the following discussion assumes there is only one 

distinct ion species of each polarity. 

There is one equation for each of the possible charge states (k,i) of an aerosol 
particle. The series of charge states is generally truncated when the probability ot 
achieving such a state becomes insignificant. Each equation has two birth terms 
corresponding to ion attachment to particles (of the same size) m adjacent charge states 
and two death terms corresponding to ion attachment to particles m the given charge 
Aerosol coagulation is generally ignored. For a polydisperse aerosol modeled as a set of 
discrete size classes, there is a set of balance equations. for each size class and the ion 
equations have a death term associated with each possible charge state and size class of an 
aerosol particle. For a polydisperse aerosol modeled as a continuum, there is one size- . 
dependent aerosol equation for each charge state and the ion-aerosol death terms m the ion 
equations become integrals over particular size. 

In addition to the time derivative equations there may be conservation equations for 
charge and/or total particles of a given size. These are generally not. all independent so 
that one or more may be eliminated. In the general case, the remaining system of 
equations must be solved simultaneously. In some cases, the.ion concentrations may be 
treated as known in which case the ion time derivative equations are no longer required 
and the set of equations for each aerosol size may be solved independently of other sizes. 
This is the case when the ion-aerosol death terms in the ion equations are small with 
respect to the ion-ion recombination/death term. The ion equations. can then be solv 
independently of the aerosol component to obtain the ion concentrations. In the case where 
ion generation is by radiation and the only significant death mechanism is by ion 
recombination then ion birth and death is always in pairs and the positive and negative ion 
concentrations must be equal (assuming there is no external field to separate the ions) 
Alternatively, the ion concentrations might be obtained by direct measurement. When 
steady-state conditions hold (time derivatives equal to zero) then the aerosol equations 
depend only on the ratio of positive and negative ion concentrations and absolute ^ 
concentrations are not required. If ion birth and death are by radiation and recombination 
alone, then the ion concentration ratio is necessarily 1 and the ion equations need not be 
solved at all. 

In most ambient stratospheric situations the total aerosol concentration is likely to 
be low enough that ion-aerosol interactions can be neglected in the ion balance equations. 

It is also likely that the steady-state solution may be used. If both conditions hold and a 
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neutralizer is not used before the DMA, then the ion concentration ratio is 1 and only the 
aerosol equations need be solved as indicated above. This may not be the case for the Hig 
Speed Research Program. It depends on the rate and extent of dilution of the aircraft 
exhaust. This also determines whether or not aerosol coagulation may be neglected m the 
charging equations. The relative magnitudes of the various birth/ death terms m the 
charging equations cannot be determined until the attachment coefficients have been 
determined. 


ATTACHMENT COEFFICIENTS 

The next step in the aerosol charge distribution calculation is to determine the 
coefficients of the birth/death terms in trie balance equations. If the ion equations are 
needed then knowledge of the ion generation rate, the ion recombination coefficient and.the 
ion-aerosol attachment coefficients is required. The particle equations require only the ion- 
aerosol attachment coefficients. 

The ion birth term is a single number, q, indicating the ion generation rate per unit 
volume. The natural ion generation rate in the stratosphere depends on pressure 
(altitude), temperature (?), solar zenith angle and perhaps other factors. In a radioactive 
neutralizer, the ion generation rate depends on pressure, temperature (?), source strengt 
and type, distance from the source and perhaps other factors. The ion # 
recombination/death term is the product of the ion recombination coefficient a, and 
positive and negative ion concentrations. The ion recombination coefficient depends on 
pressure, temperature and gas composition through the ion properties of mass and 
mobility. (Ion properties are discussed in more detail below ) There appears to be no 
theory available which can accurately predict the ion recombination coefficient from mass 
and mobility. The recombination rate must therefore be inferred from measurements [e.g. 
of generation rate and steady-state ion concentrations). If ion-aerosol interactions m the 
ion balance equations can be neglected, then the ion balance equations need not be solved 
but instead the ion concentrations may be measured. If steady-state conditions hold, then 
only the ion concentration ratio is needed. It is not clear that the temperature and/or 
pressure dependence of any of these ion variables can be accurately predic e . no , 
measurements of ion concentrations (or whatever) would have to be obtained over the 
entire range of expected sampling conditions. 

The ion-aerosol birth/death terms are each a product of an attachment coefficient, 
an ion concentration and a particle concentration. The ion-aerosol attachment coefficient, 
0-a ki ( n ) i s the flux of ions of polarity j to a particle of diameter Dp carrying * charges ot 
polarity i normalized to unit ion concentration. The ion-aerosol attachment coefficient 
depends on ion mass and mobility and on particle size, charge level, relative polarity and 
dielectric constant (see e.g. Pui, Fruin and McMurry, 1988). Besides indirect dependence 
through the above parameters, it also depends directly on temperature and mean molecular 
mass of the gas. If three-body ion trapping is important, then the ion recombination 
coefficient is required to estimate the ion-aerosol trapping distance. 

The theory for calculating an ion-aerosol attachment coefficient is rather . 

complicated. The Fuchs’ (1963) limiting sphere approach has been generally successful m 
the transition regime. The exact radius of the limiting sphere varies from author to author 
but is generally on the order of one ion mean free path greater than the particle radius. 
Outside this sphere ion movement is treated in the continuum regime accounting lor 
diffusion and electrical forces due to the image charge as well as the net charge on the 
particle Inside the limiting sphere the ion is treated in the free molecule regime, z.e. in 
free flight within the imposed electric field. The fraction of ions entering the limiting 
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SDhere which are actually captured by the particle must be calculated. This fraction and 
the exact location of the limiting sphere is where a large measure of the uncertain!} in t .e 
overall calculation arises. 

For pure coulombic attraction or repulsion, the ion is captured if and only ml i its 
traiectorv takes it to the actual particle surface. For pure image charge there is a. cnt c 

Set parameter and minimum apse (point of closest approach) r f °' h th^par ide pass 
trp iprtnrips Ion traiectories with greater impact parameters approach the particle, pass 
Ihroueh an apse and recede. Ion trajectories with lesser impact parameters have no apse 
and ultimately spiral in to the particle surface. The fraction of ions captured by the 

panicle can be calculated analytically in either case A ‘ ^further 

when there is both a net charge and an image charge. The calculation is runner 
comDlicated if a distribution of ion velocities is considered at the limiting sphere. 

Following Keefe ct al (1968), Hoppel l 1 ? 77 ) ass , un f? 

TTnwpvpr the ions are being accelerated in the electric field of the particle, ine resulting 
Velocity bias built up by an ion between collisions with neutral gas molecules persists to 
some extent after a single collision due to the generally greater mass of the ion There may 
thus be an accumulated velocity bias of the ions at the limiting sphere which is not 
accounted for M gration ove the velocity distribution to obtain the captured ion 
S In be carried out analytically for the pure 

Maxwellian velocity distribution but numerical solution is required in other cases, u g 

‘the numerical results of Keefe et al. (1968), Hoppel developed an ^^^YtsTcurTv It reduced 
appears to work satisfactorily at normal temperature and pressure. Its accuracy at reduced 

pressure and temperature is not easily extrapolated. 

In the limit of very small particles, the ion-aerosol attachment coefficient must 
approach the Ion recombination Coefficient. One problem here is that the t theory foM on- 
aerosol attachment generally treats the particle as stationary and only the >on i mmes In 
the limit the particle mass approaches the same order as the ion mass and botn must oe^ 
considered mobile. Hoppel (1977) notes that there are compensating effects in the theo } 

which tend to cancel out this error. The other problem is that • ? a ' n of® ^osTte pSarity 

recombination is three-body trapping. An ion approaching another on of o^ ©site po y 
ea ins kinetic energy as it falls into the potential well of the attractiv e electric field. 
Collisions with neutral gas molecules tend to strip it of this excess kinetic energy un i 
total^nergy* (kinetic + potential) becomes negative. At this point it no longer has enough 
kinetic energy to escape the potential well and the ions are trapped in orbit around each 
other Further collisions modify this orbit until one ion achieves escape velocity v c an 
(presumably low probability), the ions collide or an electron penetrates the potential 
barrier between the two. It is necessary then to determine the fraction of ions which, 
become trapped in orbit around a particle due to collision with one or more neutral a,r 

molecules. 

Hoppel references Natanson (1960) with regard to his treatment of threnody 
trapping. In Hoppel’s work a trapping distance is defined as that radius i from the ? ai 
(electric field) center at which the ion accumulates a given amount of excess kinetic e gj 

within a radial coordinate change equal to one mean free path. An ion /TtTflv bv the 
Fuchs’ limiting sphere is then considered captured if it is either captured direct > ) 
electric forces (impact parameter less than the critical impact parameter descri ' 

or if i t collides wi th a neutral gas molecule within the threnody trapping distance. The 
calculation of fraction captured again relies on approximations concerning the te y 
distribution and integration thereof. 

Of greater concern though is the definition of trapping distance. _ The excels energy 
from which it is defined is obtained from a similar definition of the ion-ion trapping 



distance The ion-ion trapping distance is obtained from a measured value of the ion 
recombination coefficient using an equation of Natanson relating these two parameters. 

Since the analysis only looks at the excess kinetic energy accumulated in one mean free 
path length it appears that it ignores the accumulated loss of total energy which should be 
the relevant energy level to consider. It also ignores once again the persistence of velocity 
(hence, also excess kinetic energy) of the larger mass ion. In addition, Hoppel s definition 
of trapping distance from excess kinetic energy does not account for the random 
distribution of ion impact parameters. 

Attempts to calculate three-body trapping distances directly from theory have not 
produced accurate results. Hoppel cites Brueckner (1964) as an example and uses ; his i work 
to deduce the dependence of the ion-ion trapping distance on the lon/neutral gas molecule 
mass ratio. However, Brueckner not only failed miserably m predicting the .magnitude of 
the trapping distance, but he also failed to weight his integrations over possible collision 
orientations by the probability of a collision which is proportional to the relative velocit} of 
the colliding ion and neutral molecule. 

[Consequently, from what I have read to date (which does not include Nat anson’s 
work nor any other primary reference on three-body trapping), calculations of three body 
trapping distances seem rather questionable whether they are directly from theory or rely 
on measurements of ion recombination coefficients.] 


Figure 3 in Hoppel (1977) and Figure 4 in Hoppel and Frick (1986) indicate that 
three-body trapping is only important for particles smaller than a few hundredths of a 
micron in radius at normal temperature and pressure. .An equation m Hoppel and Frick 
(1986) (p 12) gives dependence of trapping distance, d, on mean free path, A, as 
doc Af(i+5e 2 /3l:TA)^ 2 -l] (attributed to Natanson). As pressure, p, is reduced, d goes up 
far more slowly than A which is inversely proportional to p. The probability of collision 
with a neutral gas molecule within the trapping distance should go roughly as 
l-exv(-2dl\) which should therefore decrease as pressure decreases. One would reason 
then as did Hoppel (private communication), that three-body trapping decreases m 
importance at higher altitudes and reduced pressures and temperatures. Changes m d due 
to the temperature dependence should be of lesser magnitude than those due to t e 
pressure dependence. Thus, if one is willing to accept the theoretical basis of three-body 
trapping, it may be possible to ignore the effects of three-body trapping within the flight 
range of the ER-2 for particles larger than a few hundredths of a micron. 


Still of some concern is the exact location of the Fuchs’ limiting sphere. Fuchs 
himself chose the average radius attained by ions traveling a distance of one mean free path 
in a random outward direction from the surface of the particle (Hoppel and Frick, 19bo, 

Eq 16). The displacement of this sphere from the particle surface is less than one mean 
free path. (5=2 A/3 for A«aand 5=A-a/3 for A» a where sis the. displacement, A is the 
mean free path and a is the particle radius.) Hoppel chooses the limiting sphere to be one 
mean free path from the image capture radius (minimum apse), A or from the trapping 
distance, d, whichever is greater. .The sensitivity of the overall calculation of ion-aerosol 
attachment coefficient to this variation is unknown. 


ION PROPERTIES 


The last step in the aerosol charge distribution calculation is to obtain. the necessary 
ion properties. The properties needed are the mass, electrical mobility, diffusivity, mean 
thermal speed, and mean free path. Using relations cited m Hoppel and Frick (19S6) and 
Pui Fruin and McMurry (19SS), the mean thermal speed can be calculated from mass and 
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the diffusivitv from mobility at a given temperature. The diffusivity is proportional to 

mean free path, mean thermal speed and (1 + m/M) 1 ' 2 ^ ere iu^ arfes by 

the ion and neutral gas molecule, respectively. The constant of proportionality varies by 
slightly more than the square root of 2 according to whether one uses the equation of Loeb 
(Hoppel and Frick, 1986) or that of Kennard (Pui, Frum and McMurry, 1988). Using 
either of these expressions, the mean free path can be obtained from the ion mass and 
mobility and the neutral gas molecule mass. The presence of the factor containing the 
mass ratio indicates that this represents a collision-free path and not a persistence path. 

As noted above there are points in the theory where the persistence path might seem more 
appropriate A persistence path can be obtained from the Loeb or Kennard expressions by 
SBng the mass ratio factor. For m=200 amu and M=29 amu the persistence path is 
larger than the collision-free path by a factor of 2.8. 

Hoppel and Frick (1989) indicate that the calculation of attachment coefficient is 
relatively insensitive to the choice of using the Loeb or Kennard expression. Thoug 
switching to the persistence path is a somewhat greater variation the calculation may be 
relatively insensitive to this as well. Similarly, one might expect little variation due to the 
various choices of the location of the Fuchs’ limiting sphere. A more careful sensitivity 
analysis is required to verify this. 

The minimum required ion properties are thus mass and mobility for both positive 
and negative ions. Various investigators have developed empirical plots or expressions 
relating mass and mobility. See, for example, references cited m the Appendix of Hoppel 
and Frick (1990), in Adachi, Kousaka and Okuyama (1985) and m \\ uedensohler, 
Liitkemeier, Feldpausch and Helsper (1986). It is not clear how much confidence can be 
put in these expressions for the ER-2 application. 

Measurements indicate that a range of properties is to be expected for each ion 
polarity corresponding to the range of possible molecules which may be incorporated in an 
ion. Ions also undergo evolution on several time scales as mentioned m the Appendix of 
Hoppel and Frick (1990). Thus, ions produced m a neutralizer used on the EK - are likely 
to be significantly different from ambient stratospheric ions. In either case it is possible to 
treat each distinct ion species individually in the population balance equations. This o 
be quite complicated and the available data on ion properties is probably insufficient to 
support such an analysis. Most investigators have chosen to use a single set of average ion 
properties for each polarity or even the same set for both polarities though the latter option 
is not recommended. This approach appears to be generally successful but the error 
incurred is likely to increase as the range of ion properties increases. Wiedensohler and. 
Fissan (1988) attempted to account for a second negative ion species - free electrons m 
their work with ultrapure gases. It seems unlikely that there will be a significan 
concentration of free electrons in any charging situation imagined for an bit - based divla 
but this should be considered. (See also O’Hara et a/., 1989.) 

The papers mentioned in this section provide a number of references which should 
provide a good starting point for resolving the problem of determining ion properties. A 
few of these references have been already obtained and are listed under ion references. 


CONCLUSIONS 

Overall the problem of aerosol charging for a DMA operated in the stratosphere 
may range fror^ intractable to reasonable. It depends on the measurement conditions the 
charging process used, availability of ion property data, and a sensitivity analysis of the 

charging theory. 
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It is recommended that the natural charge state of the aerosol be used. Han _ 
artificial charger is used, either a radioactive bipolar charger or a corona discharge unipolar 
charger, the charging process probably cannot be considered steady-state (Hoppel and 
Frick 1990). There will be significant wall losses for the ions which will perturb the ion 
concentration ratio in a bipolar charger. There is probably far less data on a ^ Gc i y 
produced ions under stratospheric conditions than naturally occurring ones. The aerosol 
charge distribution from an artificial charger is likely to be a far stronger function of 
altitude than the naturally occurring charge distribution. 

If the naturally occurring charge distribution can be considered steady-state and if 
the only significant ion birth and death processes are radiation and recombination, then the 
aerosol^balance equations are relatively easy to solve. If three-body ion trapping can be 
neglected, then the ion-aerosol attachment coefficients can be carried out with a reasonable 
degree of confidence (subject to a sensitivity analysis). If all of the above conditions are 
met the ion generation rate and recombination coefficient are not needed. All that is 
required then are ion masses and mobilities. The question remains as to what degree of 
confidence these may be predicted from the literature. 

The conditions stated in the previous paragraph may not be satisfied when 
measuring supersonic transport exhaust. There may also be problems in measuring 
aerosols that have been in the dark for a long period. (Does the ion generation rate drop 
significantly in the dark?) If any of these conditions is not met, the complexity ot the 
calculations increases and more data is required from the literature - data that may be 
more difficult to find and less reliable than the ion masses and mobilities. In this case, our 
confidence in the calculated aerosol charge distribution is likely to drop dramatically. 



Appendix 3 




Evaluation of Diffusion Batteries 
for Stratospheric Aerosol Measurements 


Diffusion batteries have the potential to be used to measure the size distributions of 
submicron particles in the stratosphere. This report presents the results of a feasibility 
study of this application of diffusion batteries. 

A diffusion battery is device used to separate particles on the basis of their diffusional 
mobility. Particles are pulled through a “battery” of screens or tubes. Particles of a given 
size have a known probability of striking the battery matrix. This probability is higher for 
smaller particles (with higher diffusion coefficients) than for larger particles (with lower 
diffusion coefficients). By varying the flowrate and/or the number of screens or tubes, 
particles of different sizes can be selected to penetrate through the battery to reach a 
counting device (usually a condensation nuclei counter). The controlling parameter in 
diffusion battery performance is the inverse Peclet number, (j). For a tube-type diffusion 
battery, 4> is defined (e.g., Brown et ah, 1984) as 


§ = N c 7tD*L c Q \ 

where Nc is the number of tubes in the diffusion battery, D* the diffusion coefficient, Lc 
the length of the tubes, and Q the volumetric flowrate. Particle diameter enters the 
formula through the diffusion coefficient, given as 
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where k is Boltzmann’s constant, T the temperature, rj the viscosity of the fluid (air), I the 
mean free path of the fluid and D p the diameter of the particle. The experimentally- 
determined efficiency of a particle penetrating a tube diffusion battery is given by (Brown 
et ah, 1984) 

C=0.819exp(-3.650) + 0.097exp(-22.3(D) [0>0.04] and 

Q=\ -2.560 273 + 1.20 + 0.177O 473 [0<0.04]. 

Measurements are obtained from the diffusion battery by selecting a “stage”, a 
combination of flowrate, number of capillaries and capillary length. The number of 
particles penetrating each stage is then measured. These numbers (with error) and a 
knowledge of the penetration efficiency as a function of size (also with error due to the 
measurement of airflow) for each stage are used to extract the actual particle size 
distribution. The aerosol size distribution is assumed to remain constant while the number 



concentration of particles penetrating each stage of the battery is measured. If a single 
CNC is used to cycle through each stage, obtaining statistically reasonable counts may 
take several minutes to tens of minutes at stratospheric concentrations. 

A diffusion battery using commercially available glass capillary array filters has been 
modelled. The particle penetration efficiency for an ambient pressure of 50 mb at an 
instrument temperature of 273 K has been calculated as a function of particle diameter for 
each stage of the battery. These efficiencies are shown in Figure 2. The key feature of the 
penetration efficiency curves is the extremely broad nature of the response function for 
each stage of the battery. The measured response-the number concentration of particles 
penetrating each stage— is a convolution of the particle size distribution and these curves 
of particle penetration efficiency: 


Cj = J$(D)N(D)dD, (3.1) 

Di 

where Cj is the concentration of particles penetrating stage j of the diffusion battery, Q the 
penetration efficiency of particles of size D penetrating stage j, N the concentration of 
particles of size D in the aerosol, and Di and D 2 the minimum and maximum particle 
diameters detectable with the CNC. Recovering the particle size distribution from the 
measurements Cj and known response functions Q(D) involves inverting the integral to 
recover N(D). Unfortunately, for any given set of measurements Cj, an infinity of possible 
solutions exist that satisfy(3.1). Additionally, because the problem is extremely ill-posed, 
small errors in measurements can lead to very large errors in the inverted size distribution 
using traditional matrix inversion techniques. For most applications, the greatest difficulty 
in using diffusion batteries involves this inversion process. Because this process is the 
limiting factor in using diffusion battery data for stratospheric aerosol size distribution 
measurements, considerable effort was expended in examining various techniques. 

The aerosol literature is full of techniques for solving the inversion problem. Most telling 
is that no single procedure has been universally accepted, and papers extolling new 
inversion techniques for diffusion batteries continue to appear. For this study, two new 
techniques for inverting diffusion battery data have been examined. The first of these, the 
Extreme Value Estimation (EVE) procedure, has been developed by P. Paatero of the 
University of Helsinki (Paatero, 1990). This technique does not attempt to find a single, 
“correct” size distribution, based on the argument that any size distribution that is 
consistent with the integral equation is valid. Instead, the EVE technique uses matrix 
inversion techniques to finds a solution set, limited to non-zero values, that is consistent 
with the measurements and their errors. The EVE program, which is commercially 
available from Dr. Paatero, was purchased for the purpose of this study. The source code 
for the software is not available. 

The EVE program was used to invert a number of modelled size distributions. For 
modelled stratospheric measurements, the errors in the measurements due solely to 
counting statistics were sufficiently large to render the results essentially meaningless. The 



deconvoluted family of size distributions that satisfied the equations within the estimated 
errors were broad enough to prevent reasonable determination of either the mean size or 
the geometric standard deviation of a simple lognormal, monomodal distribution. 

An alternative method, the maximum entropy (MAXENT) technique, is based on the 
argument that the “best” solution of the particle size distribution is that containing the 
least “information content” in the language of information theory (Yee, 1989). The 
solution with the least information content is the simplest and smoothest of the possible 
solutions; it has no extraneous information beyond that warranted by the measurements. 
This concept can be expressed mathematically by the Shannon-Jaynes entropy equation, 


H(n) = J n(D)ln[n(D)]dD . 


Di 


where H(n) is the Shannon-Jaynes entropy, and n(D) the normalized particle size 
distribution function (i.e., the particle probalility distribution) such that J n(D)dD=l. 

Di 

The MAXENT algorithm searches for n(D) that satisfies Eqn. 3.1 and that maximizes 
H(n). This method requires no a priori, ad hoc assumptions about the expected particle 
size distribution, as do most other inversion techniques, and produces a single size 
distribution by introducing the rational and valid constraint of maximum entropy. 

The MAXENT algorithm was obtained from Dr. E. Yee of the Canadian Defence 
Research Establishment. Again, the source code was not made available. As before, a 
number of modelled diffusion battery measurements based on a known particle size 
distribution function were used as input to the inversion algorithm. Measurement errors 
were also used as input, although the penetration functions were assumed to be perfectly 
known (no allowance was made for errors in flow measurement). The MAXENT 
algorithm showed excellent potential by correctly recovering some of the input size 
distributions with very good accuracy when measurement errors were small. However, 
the algorithm was inherently unstable; tests using realistic particle size distributions with 
realistic errors did not recover any of the input distributions and instead oscillated wildly. 
These failures may be in part a product of numeric difficulties in the complex algorithm to 
solve the maximum entropy equation. Until these limitations are reduced, the algorithm is 
of little use when realistic measurement errors are introduced. 

Conclusions 

Diffusion batteries coupled with CNCs show some potential for estimating the size 
distribution of stratospheric particles. The technique is severly limited by the lack of 
robust algorithms for retrieving data with measurement errors. The MAXENT method 
shows some potential when measurement errors are small. By artificially introducing small 
errors into the MAXENT program, some qualitative information on the mean particle size 
may be obtained, although the substantial errors in the details of the retrieved size 



distribution may result. Diffusion batteries cannot now be recommended for 
measurements of the stratospheric particle size distribution. 
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